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_ ■ Abstract 
> 

Employing the closed-time path 2PI effective action (CTP 2PI EA) approach, we 
study the response of an open interacting electronic system to time-dependent ex- 
ternal electromagnetic fields. We show that the 2PI EA provides a systematic way 
of calculating the propagator and response functions of the system. Due to the 
invariance of the 2PI EA under external gauge transformations, the response func- 
tions calculated from it are such that the Ward-Takahashi hierarchy, that ensures 
current conservation beyond the expectation value level, is satisfied. These findings 
may be useful in the study of interacting electronic pumping devices, and serve to 
^ clarify the connection between current conservation (beyond the mean value level) 

and real-time nonlinear response theory. 

Key words: 2PI effective action, nonlinear driven transport, generalized 
Ward-Takahashi identities 

PACS: 05.30.-d, 05.60.Gg, 72.10.-d, 72.10.Bg, 73.23.-b 



O 



Email addresses: jperalta@df .uba. ar (J. Peralta-Ramos), 
calzetta@df.uba.ar (E. Calzetta). 



Preprint submitted to Elsevier 



13 January 2009 



1 Introduction 



In this work we determine the basic requirements that a field-theoretical ap- 
proach to open and driven systems must satisfy in order to produce current 
conserving results (in the sense of the Ward-Takahashi hierarchy developed in 
Sec. 12]), in transport calculations going beyond linear response. Another aim 
we have in mind is to clarify the close relation that exists between current 
conservation and response theory, especially in the nonlinear regime. In order 
to analyze these issues, we combine the so-called external gauge invariance 
method with the closed-time path two-particle irreducible coarse-grained ef- 
fective action (CTP 2PI CGEA), which has proved a valuable tool in the 
study of strongly interacting quantum open systems both in and out of equi- 
librium [1,2]. Although our analysis is quite general and can be applied to a 
variety of condensed matter systems, we have in mind a system of strongly 
interacting electrons subjected to external driving fields and in contact with 
two ideal reservoirs of noninteracting electrons. This could be a picture of the 
so-called interacting electron pumping devices, which are currently attracting 
much interest both experimentally and theoretically [3]. 



2 Ward-Takahashi hierarchy 



In many-particle systems, n-point vertex functions represent generalized cur- 
rents which satisfy the hierarchy of Ward-Takahashi identities [4]. This hi- 
erarchy is satisfied to all orders in the exact theory, warranting local gauge 
invariance and the conservation of the associated charges. In its most gen- 
eral form, the generalized continuity equation (WTH) can be written as (we 
use the simplified notation 1 = (£i,ri) and employ Schwinger-Keldysh non- 
equilibrium formalism [1,5,6]): 

^Af n) (l...n, l'...n'-z) = 

i n e{-5(z - n)G n {\ ... n, 1' ... (n - l')z) + ... 

+ {-l) n 5{z - V)G n {\ ...n,2' ...z) + ... (1) 

+ 5(z - n)G n (l . . . (n - l)z, l'...n')-... 

-(-l) n 8(z-l)G n (2...z,l'...n')} 

where 

A» n) (l2...n,l'2'...n';z) = 

< 1^(^(1)^(2) . . . ^(n)V>V) . . . ^{2')^{l') > 
is the (n + l)-point vertex function with current insertion at z = (t z , r z ), 
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f(z) = -e lim D IJ '(z, z')ip\z')ip(z) 



D\z,z') = (2 i )- 1 (Vl-Vl, 

-ie[A\z) + A\z')}) n 
D°(z,z r ) = 1 



1,2,3 
= 



(3) 



with the external fields denoted by A^, and G n are real-time propagators 
defined as usual 



The classical continuity equation and the usual WT identity correspond to 
the cases n=0 and n—1 of Eq. (CQ), respectively. 

As the WTH shows, in the exact theory particle number is strongly conserved, 
not only in the mean. This means that for a many-particle system in the 
presence of an external field, the current is conserved to all orders in the 
external perturbation. 

It is impossible, in general, to obtain the propagators of an interacting field 
theory exactly, and some approximations must be made. In particular, for a 
strong coupled theory re-summation concepts are usually needed [1,7], and 
warranting conservation laws then becomes a nontrivial issue. A systematic 
way of generating conserving approximations (at the classical level, i.e. n=0) 
was given by Baym [8] and corresponds to his well-known ^-derivable scheme. 
The solutions obtained from truncating the $ functional are such that the 
expectation values of the respective Noether currents are conserved [1,8]. We 
emphasize that this situation corresponds to the case n=0 of the WTH. There- 
fore, the conservation of generalized currents encoded in the WTH is not au- 
tomatically warranted in this approach. 

Although in some cases the conservation of current at the expectation value 
level may be sufficient, it is not suitable for nonlinear response studies. The ex- 
ternal gauge method [9] is a possible way of overcoming this problem present in 
^-derivable approaches. It provides a systematic way of generating consistent 
Scwhinger-Dyson (SD) and Bethe-Salpeter (BS) equations that automatically 
satisfy WT identities. Most importantly in the context of transport theories of 
strongly interacting systems, the derivation of the SDE and the BSEs can be 
done to any order in the external field coupled to the system. Therefore, the 
latter is especially suited for the study of response theory beyond first order 
and its relation to current conservation. 

1 Another method is the re-summation technique due to Hees and Knoll [10], but 
we shall not discuss it here. 



G n (l...n, l'...n') = 

r (n) < T c i/j(1) . . . ^(n)^(n) . . . 4)\l') > 



(4) 
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3 CTP 2PI effective action 



The system we are interested in consists of a central region with interacting 
electrons described by fields ijp\tp) and coupled to an external field. This 
central region is connected to two ideal reservoirs of noninteracting electrons 
described by fields ((/>' a , <f> a ) , with a = (L,R) denoting left or right reservoir. 
The reservoirs are assumed to remain in equilibrium at all times. The CTP 
classical action [1] of the system is S[ip, ip, 0, (f)\ — + + S c , with: 

= c AB ^ A B^ B <t> B (5) 
S c = c AB i) A T AB (j) B + h. c. 

where 

being Uabcd the completely antisymmetrized bare interaction local vertex. 
Tab is a local coupling parameter between the central region and the reservoirs. 

We are using a DeWitt notation [1] with A = (x,a), x = (t x ,r x ,a) and 
a = (1,2) [or (+,—)] being CTP indices indicating the branch within the 
closed-time contour. For the fields describing electrons inside the reservoirs, 
an additional index a must be included in the CTP indices (A, B), but for 
simplicity we leave it implicit. Repeated indices are assumed to be integrated 
or summed. c a b is a CTP metric c a b = diag(l, —1), while cab = c a b5(xA,XB)- 
Cab and B A b are the free CTP propagators corresponding to and S^. They 
satisfy equations of motion governed by single-particle Hamiltonians 

ho(l) = J[-iVi + eA(l)] 2 + eA (l) 

\ (7) 
VI) = - gV? , 

respectively. By performing the double Legendre transform on the generat- 
ing functional of connected propagators, and then using the background field 
method we can write the CTP 2PI CGEA of the system as [1] 

F[G, A] = iTr InG - iD AB G AB + T 2 [G] , (8) 

where Dj B = C AB -HX^ab, being E^ab the self-energy due to the reservoirs. 
T 2 [G] encodes all quantum corrections and consists of vacuum 2PI closed di- 
agrams with full propagators in internal lines and vertices corresponding to 
a theory with shifted classical action S[ip + ■?/>, ^ + ip] (neglecting constant 
and linear terms), where denote fluctuations. The diagrams are vac- 

uum because the mean value of the fluctuation fields is zero by construction. 
Note that because of the vanishing of the mean fields, the shifted classical 
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action vanishes when evaluated at (ip, ip). The case of noninteracting electrons 
corresponds in this scheme to T 2 = 0. 

The Schwinger-Dyson equations for the CTP propagators follow directly from 

r[G,A] 

sr 1 

J7T- = —K AB (9) 
oG AB 2 

where the physical case corresponds to vanishing external sources K — 0. 
Defining the self-energy 

Z AB [G] = -25T 2 [G,A]/5G AB (10) 
we can rewrite the SDE in the usual way 

G AB = D AB + i^ AB [G] , (11) 
being T, A b one-particle irreducible by construction. 

We note that in the 2PI EA approach, mean fields and two-point propagators 
are treated on the same footing, whereas higher order propagators are ob- 
tained from them [1]. Although this represents an approximation to the exact 
dynamics of the quantum system, we will show that it is consistent with the 
WTH. 



4 WT hierarchy from external gauge invariance 



Under a local transformation U(l) — expieip(l), the external field and the full 
propagator transform as [2, 11] (we omit the CTP indices for the moment) 

A, -> 4 = UA.U- 1 - i{d»U)U~ l 

G(l,2) ^ G'(l,2) = U(1)G(1,2)U- 1 (2) . 1 ' 

The crucial observation that allows us to relate the CTP 2PI CGEA of the 
open system to nonlinear transport through it is that T[G, A] is invariant under 
a gauge transformation of the external field A M . Following Bando, Harada and 
Kugo (Ref. [9]), we will call this external gauge invariance (EGI) of the 2PI 
EA. It is rather straightforward to prove that, if we retain all terms in the 
loop expansion of r 2 , then the 2PI CGEA is EGI [1,2,11]. This also holds if 
we truncate the loop expansion of T 2 at arbitrary order [12]. 

We will now make a connection between EGI and nonlinear response. For 
notational simplicity, CTP indices are omitted in the following. The solution 
G[A\ to the SDE can be expanded in powers of the external field A^ (see 
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Refs. [6,13]): 



G[A] (X, Y) = G[0] (X, Y) + iA ll B$+ 



^A^A^Tl^ + — A^A u A p Tl^ 



i ? '3 
11 1-1-111/ " 



= J2-fdl...dnA, 1 (l)...A, n (n) 



(13) 



n=0 



U^(l,...,n;X,Y) 



where in the second line we have made explicit the internal (integration) 
(1,. . . ,n) and the external (X,Y) variables. The "response" functions Hr n +2) 
encode the variation of the full propagator with the external field. It is worth 
remarking that the use of the closed-time path method automatically ensures 
that the response functions are causal [1,6]. 

Returning now to Eq. (TT3"]) . the response functions are given by 



and correspond to (n + 2)— point functions with n current vertices inserted at 
locations (1, . . . ,n) where interactions between the current and the external 
classical field take place. The functions il( n+2 ) are obtained from the SDE by 
functional differentiation with respect to A^ (and then setting A = 0). This 
results in the BSEs for the response functions. 

The combination of the SDE and the BSEs completely determine the full 
propagator and the response functions II( n+ 2). The SDE is obtained from the 
2PI CGEA, while the BSEs are obtained from the SDE by differentiation 
with respect to the external field. The important point is that, because the 
2PI CGEA is invariant under external gauge transformations, both the full 
propagator and response functions obtained this way are external gauge co- 
variant. As we will show below, the EGI property of the 2PI CGEA implies 
that G and ri( n+ 2)> as obtained from T[G, A], satisfy the WT hierarchy. This 
provides the required link between current conservation in nonlinear response 
and the external gauge invariance of the 2PI CGEA, and also a powerful and 
systematic way of studying nonlinear response in strongly interacting systems 
coupled to ideal reservoirs. 

To see the connection between EGI and the WT hierarchy, recall that external 
gauge invariance of the effective action means G[A'} = UGIA}!/' 1 . Inserting 
the expansion in powers of the external field given in Eq. (1131) into both sides 




(14) 
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of this identity we get 

g[o] • • - i;,,i;,n 



i 2 



(15) 

UGlOp- 1 + iA^UU^U- 1 + -A^UUtU- 1 + ... 



In particular, for an infinitesimal external gauge transformation U(X) wl + 
ieip(X) the transformed external field is 

= A,(X) + d^{X) , (16) 

so Eq. ( fl5l) becomes 



(17) 



E -f / [dn] m^M) + d„<p(i)}V$$ = 

ra=0 ■ i=l 

£ S / H ft 4(^){nK + *(x)ng;-y- 

« — n 'i* ^ „• — 1 



n=0 ,L - ° i=l 

where we have suppressed the arguments of nj^'^" and defined [dn] = dl . . . dn 
for brevity. Comparing terms of the same order in A^ on both sides of this 
expression we get 

n 2 (x,r) + * J di[9>(i)M(x,F;i) = 

U 2 (X,Y)+ie[ip(X)U 2 (X,Y)-U 2 (X,Y)^(Y)\ 

for the zeroth order term. Higher order terms have a similar structure but 
can become quite involved [12]. The main point to emphasize is that the EGI 
of the 2PI CGEA implies relationships among the response functions, shown 
explicitly in Eq. ffTTl) . 

Assuming that (p vanishes at infinity, we can integrate by parts the second 
term in Eq. (|18|) (zeroth order in A^) to obtain 

l JdMl)dlmX,Y;l) = 

- ie[<p(x)n 2 {x, y) - n 2 (x, y) v {y)\ 



which implies 

dl < T c f{z)^{X)i,{Y) >= 



e < T c ^{X)ijj{Y) > [5{Y - z) - 6(X - z)] 

This is precisely the identity corresponding to n = 1 in the WT hierarchy 
given by Eq. ([T]). We see that, even at zeroth order in the external field, the 
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relation between the three-point vertex and the two-point function, Eq. ( |20|) . 
is satisfied due to the EGI of the 2PI EA. 

It is clear that this procedure could be continued to higher order terms in 
A^, thus generating higher order WT identities. We note that the hierarchy 
obtained for the response functions II( n +2) is completely equivalent to that 
involving A( n ), given in Eq. ([1]), as expected since, ultimately, they both enforce 
current conservation [12]. 



5 Conclusions 

We have determined the basic requirements that an approximation to a non- 
equilibrium many-body problem in an open and driven fermionic system must 
satisfy in order to achieve current conservation beyond the expectation value 
level. One of the most important results of this work is the close relation 
found between nonlinear response theory and the Ward-Takahashi hierarchy, 
necessary for current conservation. This connection was clearly displayed by 
using an approximation scheme based on the CTP 2PI CGEA, that provides 
equations for the propagator and response functions consistent with the WTH. 
We emphasize that, although being an approximation to the true dynamics 
of the system, the 2PI EA description is consistent with the WTH, even in 
the nonlinear regime. Our results may be of use in the theoretical study of 
quantum transport through interacting electronic pumping devices, which are 
nowadays receiving much attention. 

We thank Liliana Arrachea, Carlos Naon, Alfredo Levy Yeyati and Sergio Ulloa 
for useful discussions. This work has been supported in part by ANPCyT, 
CONICET and UBA (Argentina). 
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